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I — ' Lyman break galaxies (LBGs) at z ^ 3 exhibit number densities and clustering similar to 

local i* galaxies; however, their star formation rates (SFRs) are much higher. We explore 
the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an 
explanation for these key properties. The relative paucity of starburst galaxies at low redshift can 

f^ ', be attributed to a much higher collision rate in the past. We use high-resolution cosmological A^- 

body simulations and a hierarchical halo finder to estimate the galaxy collision rate as a function 
of time in the popular ACDM cosmological model. We find that bright coUisional starbursts are 

Q^ I frequent enough to account for most of the high-z (2.5 — 4.5) LBGs. Moreover, many of the 

0^ ' objects are of relatively small mass, but they cluster about large-mass halos. They therefore 

exhibit strong clustering, similar to that observed and stronger than that of the relevant massive 
Ph! halos. 
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"ti I Subject headings: cosmologyitheory - dark matter ~ large-scale structure of universe, galaxies: 

C^ ' interactions - starburst, methods: numerical 
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Vh ' 1- Introduction 

C^ , 

Data from the Hubble Deep Field (reviewed in Dickinson 1998) and ground-based telescopes (Steidel 
et al. 1996a, Steidel et al. 1996b, Lowenthal et al. 1997, Steidel et al. 1998) have revealed a population of 
galaxies at z ~ 2.5 — 4.5. They were found by multicolor photometry exploiting the characteristic spectral 
attenuation shortward of the Lyman limit, so these galaxies are referred to as "Lyman break galaxies" . The 
LBGs are observed to be abundant and highly clustered, with comoving number densities and clustering 
properties at z ~ 3 comparable to those of present-day bright (> L*) field galaxies (Giavalisco et al. 
1998, Adelberger et al. 1998). They are forming stars at a high rate, comparable to local "starburst" 
galaxies, and are much smaller than similarly bright galaxies nearby. 

There are competing views regarding the nature of LBGs. In one view (Steidel et al. 1996b, Adelberger 
et al. 1998, Baugh et al. 1998, Governato et al. 1998 ,Mo, Mao, & White 1999), most LBGs are large galaxies 
quiescently forming stars at the bottom of the potential wells of massive dark matter halos. We refer to 
this idea as the "central quiescent" scenario. An alternate view (Lowenthal et al. 1997, Somerville, Primack, 
& Faber 1998) maintains that LBGs are mainly galactic starbursts triggered by collisions between small, 
gas-rich galaxies. We refer to this idea as the "coUisional starburst" scenario. 



Figure 1 can also be retrieved from http://physics.ucsc.edu/users/tsafrir/figl.html (USA) or 
http://www.fiz.huji.ac.i1/~ tsafrir/figl.html (Europe) 

Fig. 1. — Part of the simulation volume a,t z = 2.9. Spheres are dark-matter halos, color-coded by mass [the 
color bar is labeled with log(Mhaio/^~^Af0)] and sized in proportion to their virial radii. Stars indicate the 
locations of all halo collisions occurring at 3.9 > z > 2.9, an interval of ~ 600 Myr. Since each burst lasts a 
relatively short time TvIs ^ 100 Myr, only ^ 1/6 of the bursts would be seen at any one time. Most collisions 
occur in or near relatively massive halos, resulting in strong large-scale clustering. 

Within the central quiescent scenario, there is roughly a one-to-one relationship between the LBGs and 
massive halos (Stcidcl ct al. 1996b, Adelberger et al. 1998). This helps explain the strong clustering of 
observed LBGs, since in CDM (Blumenthal et al. 1984) theories of hierarchical structure formation, massive 
objects are more clustered than low-mass objects. It has been shown that the clustering properties of LBGs 
at 2 ~ 3 can be reproduced within various CDM models if the LBGs are associated with the most massive 
collapsed dark-matter (DM) halos at that epoch (Jing & Suto 1998, Wechsler et al. 1998, Adelberger et al. 
1998). More detailed modeling appears generally consistent with the central quiescent framework (Baugh 
et al. 1998, Governato et al. 1998). 

Until the simulations reported here, there have been no predictions for clustering properties and only 
rough estimates of number densities of LBGs within the coUisional starburst picture. Somcrville, Primack, 
& Faber (1999) used a semi-analytic treatment in order to compare the properties of individual galaxies in 
the two scenarios, and argued that the high star formation rates, small emission line-widths (Pettini et al. 
1998) (~ 70 km/s), young ages (Sawicki & Yee 1998) and high star formation surface densities (Lowenthal 
et al. 1997, Heckman et al. 1998) of LBGs are more easily explained within the coUisional starburst model. 
However, there were many unanswered questions since only simple approximations (Makino & Hut 1997) 
were used to estimate the merger rate of sub-halos at high redshift z ^ 3 and there was no way to calculate 
the spatial distribution of coUisional starbursts. 

In order to establish whether coUisional starbursts are a plausible origin for LBGs, we now ask whether 
the collisions in hierarchical scenarios can match the observed number density and clustering properties. 
These questions are addressed in this letter using high resolution A^-body simulations. We first present the 
simulations and accompanying halo finder, and briefly explain our prescription for identifying collisions and 
estimating the resulting luminosity. The number density of the observed LBGs is then compared to the 
simulated number density as a function of redshift, and the correlation function for the simulated LBGs is 
presented. 



2. A-body simulations and halo/collision finders 

Only recently have cosmological A-body simulations reached the stage where halo sub-structure can 
be resolved (e.g., Klypin et al. 1999). Our simulations make use of the ART code (Kravtsov, Klypin, 
& Khokhlov 1997) which utilizes an adaptive grid to obtain the unprecedented resolution necessary for 
identifying collisions between well-resolved galactic halos (or sub-halos) in a cosmological volume. The 
simulations followed the evolution of the DM in the popular ACDM model in which the present mass density 
of the Universe is flm ~ 0.3, and a cosmological constant fl\ = 1 — flm = 0.7 makes the Universe flat. The 
Hubble constant is h ~ 0.7 (i.e., Hq = 70 kms~^ Mpc^^), the primordial fluctuations are scale invariant, and 
the mass fluctuation amplitude is as = 1-0. Our main results are based on a simulation run down to z = 1.7 



with particle mass nip = 1.3 x IO^H^^Mq in a 30 h^^Mpc box, but wc have also used another ART simulation 
with the same number of particles (256'^) in a 60 h^^Mpc box run to z = 0. We locate DM halos using a 
bound-density-maximum halo finder (Klypin et al. 1999), now extended to cope with halo interactions and 
sub-structure (Bullock 1999). Halos are modeled by the density profile (Navarro, Frenk, & White 1996) 
/Onfw(^) = Ps/[{f/Rs)i^ + f/Rs)'^]i which provides a characteristic radius (i?s), virial radius (i?vir) and mass 
(Mvir), and the associated fitting errors. The profile of a sub-halo may be truncated short of its i?vir- We 
treat only halos with Mvir > 7 x IO^H'^Mq (> 50 particles). This is sufficient to resolve collisional starburst 
LBGs according to the luminosities assigned by our prescription (see below). Halo collisions are identified 
using pairs of stored simulation outputs at redshifts Z2 > zi. For each halo at zi, we search for sets of 
particles that originated in different halos at Z2 . If the centers of two such sets overlap within their Rs radii 
at zi, a colhsion is declared (Kolatt et al. 1999). Figure 1 shows the locations of such collisions. 



3. Prom collisions to bursts 

We assume that each small DM halo at z ~ 3 contains a gas-rich galaxy at its center and that each 
collision results in a starburst. Mergers are included not only between isolated halos, but also between 
"sub-halos" — halos that reside within the virial radius of larger halos. Note that the analytic predictions 
for merger rates obtained by Press-Schechter approximations (Lacey & Cole 1993) are not sufficient here 
because they miss sub-halos of isolated systems, they are limited in predicting the progenitor mass spectrum 
(Somerville & Kolatt 1999), and they ignore spatial correlations. 

The top panel in Figure 2 shows the time-evolution of the number density of identified collisions as a 
function of look-back time, assuming a duration of visibility TvIs — 100 Myr. About one half of the collisions 
at z ^ 4 involve sub-halos; such collisions would have been missed without these very high-resolution 
N-body simulations, which significantly supersede the cruder treatment of collisions in the semi-analytic 
investigation (Somerville et al. 1999b). Unbound collisions (those in which the two halos are not subsequently 
bound) are not accounted for at all in the semi-analytic models. Also, the dynamical friction time scale used 
in SAMs tends to over-estimate the host - subhalo collision time scale (Kolatt et al. 1999), the subhalo 
- subhalo collision approximation (Makino & Hut 1997) is somewhat simplistic, and the progenitor mass 
spectrum is uncertain (Somerville & Kolatt 1999, Somerville et al. 1999a). The collision rate per physical 
volume (not shown) declines oc (1 + z)^ for z ^ 2.5, in general agreement with theory (Lacey & Cole 1993) 
and observations. This was determined for halos > 10^^/i^^Mq using our simulation in a larger volume with 
lower mass resolution. 

To assess how many collisions should actually be observable, luminosities are assigned as follows. We 
assume that before a collision each galaxy has a cold-gas reservoir rUg = fgfhiTihaio, where /b is the fraction 
of mass in baryons (/b = fib^m^) and /g is the fraction of baryons in cold gas. We assume fib = O.OlSh^^ 
and /g = 0.3. Based on simulations including gas dynamics and star formation (Mihos & Hernquist 1994a, 
1994b, and our new simulations at z ^ 3 with /g ~ 0.3 using an updated version of the same code), we 
divide the collisions into major {m2/mi > 0.25) and minor collisions, and assume that during a burst of 
duration Tburst = 50 Myr, 75% and 50% of the gas is converted into stars respectively. Gas depletion due to 
multiple collisions is ignored because we find that at most 4% of the matter in colliding halos at z > 2 has 
participated in a previous encounter. We estimate the apparent magnitude of collisional starbursts A^ab 
in the band equivalent to 1600 A rest-frame, and TvIs, using Bruzual-Charlot (GISSEL98) stellar-population 
synthesis models (assuming solar metallicity and a Salpeter initial mass function). The upper heavy line 
in the bottom panel of Figure 2 shows the time-evolution of the number density of observable LBGs with 
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Fig. 2. — Comoving number density of collisions and observed LBGs as a function of look-back time. The 
solid line in the top panel refers to halo collisions in an TV-body simulation of ACDM. Here collisions between 
halo with masses > 7x lO^ft-~^M0 are assumed to be visible for a constant duration Tvis — 100 Myr. Errors are 
estimated from the collision-identification procedure. The dashed curve is the partial contribution (~ 50%) 
from collisions involving at least one sub-halo and the dotted line is the contribution of unbound collisions. 
The bottom panel compares the number of LBGs actually observed with predictions. The thick solid curve 
represents collisional starbursts brighter than apparent magnitude A^ab = 25.5, corresponding to rest-frame 
1600 A. Observational estimates for LBGs (TWab < 25.5) from the Hubble Deep Field (Pozzetti et al. 1998) 
and newer ground-based observations with a much larger survey volume (Adelberger et al. 1998, Steidel et 
al. 1998) are depicted by the open circles and solid diamonds respectively. The upward arrows on the data 
are an approximate correction for inferred dust absorption in the LBGs (see text). The thin solid line is 
estimated for the central quiescent scenario (see text). 



A4ab < 25.5, for the collisional starburst model. At z ~ 3, only burst events involving halos > 8 x IO^Mq 
contribute to the population of M-ks, < 25.5 galaxies (> 4 x lO^M© at z '^ 2). Because of compensating 
effects, varying both Tburst and star formation efficiency by factors of 2 either way results in changes < 30% 
in observable number density. The observed number densities of LBGs brighter than A^ab — 25.5 are shown 
for comparison, calculated from the latest data (Pozzetti et al. 1998, Adelberger et al. 1998, Steidel et al. 
1998) assuming the simulated cosmology. The predicted number densities are somewhat larger than those 
observed, and thus allow for dust extinction, which we have not included. The arrows on the data points 
result from assuming a (conservative) factor of three in dust extinction, coupled with the z ^ 3 and z ~ 4 



luminosity functions estimated by recent ground observations (Steidel et al. 1998). 

For comparison, we show the predictions of a central quiescent model. We assume that every sufficiently 
massive halo hosts one LBG whose luminosity is tightly correlated with the halo mass, and we obtain 
an effective (constant) mass-to-light ratio (M/L) by adjusting the halo mass threshold at z = 3 (M ~ 
8 X lO^^h^^MQ) to reproduce the observed abundance of LBGs (Steidel et al. 1998). We then predict the 
density evolution assuming that M/L is constant with redshift (note that a larger mass threshold is required 
at higher z for a fixed apparent magnitude limit). This type of model generically predicts a steep fall-off in 
density towards higher z, in contrast with the coUisional scenario, and in apparent disagreement with the 
newest data (Steidel et al. 1998) (solid diamonds). Redshift evolution is thus a key discriminant between 
the scenarios. The predictions of the simple model shown here are similar to those of more detailed central- 
quiescent models (Baugh et al. 1998, Somerville et al. 1999b), but large uncertainties remain. Very different 
results can be obtained if the efficiency of star formation varies with redshift (Somerville et al. 1999b). 
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Fig. 3. — Two-point autocorrelation functions. The upper curve is for collisions between DM halos of 
Af > 7 X 10^h~^MQ that occurred in the redshift interval 3.9 > z > 2.9. The error bars are combined 
Poisson and model-fit (i?s) errors. The middle curve refers to halos at z = 2.9 with Al > lO^'^h^^AlQ, and 
the lower curve is for the underlying dark matter. The coUisional starburst scenario thus predicts that LBGs 
should be strongly clustered, more than the halos themselves. 



4. Clustering properties 

A key observed statistical property of the LBGs is the strong clustering they manifest. In the central 
quiescent model the number density of LBGs sets their mass scale. Analytic approximations (Mo et al. 
1999) can then be used in order to derive their clustering properties. Here, since the collisions are selected 
by their dynamics, one must calculate clustering properties directly from the simulations. 

Figure 3 depicts the correlation function of the collisions. In the range 1 — 5 /i^^Mpc, it can be ap- 
proximated by a power law, £,c{r) ~ {r/ro}^^, with tq ~ 5 h^^Mpc and 7 ~ 2.6. Shown for comparison 
is ^h{r) ~ (r/3.5/i^^Mpc)^^'^ for halos > IO^^H^^Mq-, these are what the central quiescent scenario would 
identify with LBGs. Given current uncertainties, both correlation functions are consistent with the param- 
eters derived from observations (Adelberger et al. 1998, Giavalisco et al. 1998) for the simulated cosmology: 
tq ^ 6 h^^Mpc and 7 ~ 2. Also shown for reference is Cm('') of the underlying dark matter. The relative 
biases at 5 h^^Mpc are (Cc/Cm)^/^ ^ 3 and (Ch/^m)^/^ ^ 2. We find that -- 30% of the collisions at z - 3 
occur within halos of mass ^ lO^^'^/i~^M0, and that most of the rest occur in dense environments near 
such halos. This is consistent with the fact that the collisions are highly correlated and is confirmed by the 
similarity between the auto-correlation of collisions and the cross-correlation of collisions and halos of mass 
> lO^^°/i^^M0 (at 0.5 < r < 5 /i^^Mpc). The combination of gravitational collapse and a collision rate 
proportional to the square of the halo number density can explain why collisions occur mainly near the most 
massive halos. A visual confirmation of this trend can be gathered from Figure 1. 



5. Conclusions 

The results presented here provide the first quantitative results on the clustering of colliding halos at high 
redshift and a much more accurate measure of their number density than earlier semi-analytic calculations 
(Somerville et al. 1999), which these new results generally confirm. In particular, they show that starbursts 
associated with collisions of relatively low mass halos are consistent with the observed number density and 
clustering of bright LBGs at z ^ 2.5 — 4.5, at least in the ACDM cosmology; other popular cosmologies 
should be qualitatively similar (Somerville et al. 1999b, Wechsler et al. 1998). Finally, they predict that 
a key test between models is the number of LBGs vs. redshift, which falls off much faster in the central 
quiescent scenario than for coUisional starbursts. In addition to further tests, such as the luminosity function 
and virial mass measurements, this should ultimately distinguish between these scenarios. 

The simulations were performed at NRL and NCSA. This work was supported by grants from NASA and 
NSF at UCSC and NMSU, and by Israel Science Foundation and US-Israel Binational Science Foundation 
grants. P.J. acknowledges a Sweden-America Foundation Graduate Fellowship, and J.R.P. acknowledges 
a Forchheimer Visiting Professorship at Hebrew University. We thank Stephan Chariot for letting us use 
GISSEL98. 



REFERENCES 

Adelberger, K., Steidel, C, Giavalisco, M., Dickinson, M., Pettini, M., & Kellogg, M. 1998, ApJ, 505, 18 
Baugh, C, Cole, S., Frenk, C, & Lacey, C. 1998, ApJ, 498, 501 
Blumenthal, G., Faber, S., Primack, J., & Rees, M. 1984, Nat, 311, 517 



Bullock, J. 1999, Ph.D. thesis, University of California, Santa Cruz 

Dickinson, M. 1998, in S. F. M. Livio, P. Madau (eds.), The Hubble Deep Field, Cambridge: Cambridge 
University Press 

Giavalisco, M., Steidel, C, Adelberger, K., Dickinson, M., Pcttini, M., & Kellogg, M. 1998, ApJ, 503, 503 

Governato, F., Baugh, C, Cole, S., Lacey, C, Quinn, T., & Stadel, J. 1998, Nature, 392, 359 

Heckman, T., Robert, C, Leitherer, C, Garnett, D., & van der Rydt, F. 1998, ApJ, 503, 646 

Jing, Y. & Suto, Y. 1998, ApJ, 494, L5 

Klypin, A., Gottlober, S., Kravtsov, A., & Khokhlov, A.M. 1999, ApJ, in press 

Kolatt, T., Bullock, J., Sigad, Y., Primack, J., Dekel, A., Kravtsov, A., & Klypin, A. 1999, in preparation 

Kravtsov, A. V., Klypin, A. A., & Khokhlov, A. M. 1997, ApJS, 111, 73 

Lacey, C. & Cole, S. 1993, MNRAS, 262, 627 

Lowenthal, J., Koo, D., Guzman, R., Gallego, J., Phillips, A., Faber, S., Vogt, N., lUingworth, G., & 
Gronwall, C. 1997, ApJ, 481, 673 

Makino, J. & Hut, P. 1997, ApJ, 481, 83 

Mihos, J. & Hernquist, L. 1994a, ApJ, 425, L13 

Mihos, J. & Hernquist, L. 1994b, ApJ, 431, L9 

Mo, H., Mao, S., & White, S. 1999, MNRAS, 304, 175 

Navarro, J., Frenk, C, & White, S. 1996, ApJ, 462, 563 

Pettini, M., Kellogg, M., Steidel, C, Dickinson, M., Adelberger, K., & GiavaHsco, M. 1998, ApJ, 508, 508 

Pozzetti, L., Madau, P., Zamorani, G., Ferguson, H. C, & Bruzual, A. G. 1998, MNRAS, 298, 1133 

Sawicki, M. & Yee, H. 1998, AJ, 115, 1329 

Somerville, R. & Kolatt, T. 1999, MNRAS, in press 

Somerville, R., Lemson, G., Kolatt, T., & Dekel, A. 1999a, MNRAS, accepted 

Somerville, R., Primack, J., & Faber, S. 1999b, MNRAS, in press 

Steidel, C, Adelberger, K., Giavalisco, M., Dickinson, M., Pettini, & M. a. 1998, ApJ, submitted, astro- 
ph/9811399 

Steidel, C, Giavalisco, M., Dickinson, M., & Adelberger, K. 1996a, AJ, 112, 352 

Steidel, C, Giavalisco, M., Pettini, M., Dickinson, M., & Adelberger, K. 1996b, ApJ, 462, L17 

Wechsler, R., Gross, M., Primack, J., Blumenthal, G., & Dekel, A. 1998, ApJ, 506, 19 



This preprint was prepared with the AAS lAT^rjX macros v4.0. 



This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/9906104vl 



